Current-induced domain wall motion with adiabatic spin torque only in cylindrical 

nanowires 



O 



(N 



CO 



Z. Z. Surfl and J. Schliemann 

Institute for Theoretical Physics, University of Regensburg, D-93040 Regensburg, Germany 

P. YarQ and X. R. Wang 

Physics Department, The Hong Kong University of Science and Technology, Clear Water Bay, Hong Kong SAR, China 

(Dated: March 23, 2010) 

We investigate current-driven domain wall (DW) propagation in magnetic nanowires in the frame- 
work of the modified Landau-Lifshitz-Gilbert equation with both adiabatic and nonadiabatic spin 
torque (NAST) terms. Contrary to the common opinion that NAST is indispensable for DW 
motion [l], @], we point out that adiabatic spin torque (AST) only is enough for current-driven 
DW motion in a cylindrical (uniaxial) nanowire. Apart from a discussion of the rigid DW motion 
from the energy and angular momentum viewpoint, we also propose an experimental scheme to 
measure the spin current polarization by combining both field and current driven DW motion in a 
flat (biaxial) wire. 

PACS numbers: 75.60.Jk, 75.75.-c, 85.75.-d 
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The domain wall (DW) motion in magnetic nanowires 
has recently attracted much attention in the field of 
nanomagnetismjlldue to the enormous potential indus- 
trial applications |3j, such as memory bits and logic de- 
vices. Besides the field-induced DW motion[J-|7[, the 
current-driven magnetization reversal in both magnetic 
multilayers and nanowires through the spin torque (ST) 
transfer mechanism[8J has been massively studied, chiefly 
under the aspect of low power consumption and local- 
ity of electric currents. Moreover, a lar ge n umber of 
theoretical and experimental studies [{j. Ildl| were de- 
voted to the current-driven DW motion. Here one distin- 
guishes between the adiabatic spin torque (AST) which 
origins from the polarization of the itinerant electrons 
adiabatically following the magnetization direction, and 
the nonadiabatic spin torque (NAST, often referred as a 
/3-term) due to the mismatch of the current polarization 
and magnetic moments. The common current opinion in 
the literature is that the latter mechanism is indispens- 
able for current driven DW motion [l|, @j- 

Experimental studies so far have actually been done on 
flat nanowires exhibiting a biaxial magnetic anisotropy 
rather than cylindrical wires realizing the uniaxial case. 
This is mainly due to difficulties to produce cylin- 
drical metallic wires by conventional lithography al- 
though there are techniques available to grow them via 
electrodeposition[ll| . In this letter, starting from an an- 
alytical analysis of the modified Landau-Lifshitz-Gilbert 
(LLG) equation, we show that AST only is enough to 
drive a sustained DW motion in a cylindrical wire. Our 
result relies on the fact that the magnetic moments in a 
cylindrical wire precess around the easy axis, while they 
rotate in a plane in the biaxial case[5(. Furthermore, in 
the uniaxial case no Walker breakdown occurs for suffi- 
ciently large currents (or fields), as also shown recently 
by micromagnetic simulations |12j. Finally, we also pro- 



pose a scheme to measure the spin current polarization P 
by combining both field and current driven DW motion. 

A magnetic nanowire can be described as an effectively 
one-dimensional (ID) continuum of magnetic moments 
along the wire axis direction. Magnetic domains are 
formed due to the competition between the anisotropic 
magnetic energy and the exchange interaction among 
adjacent magnetic moments. Without loss of general- 
ity, we shall consider a head-to-head DW structure, as- 
suming the easy axis being along the wire (z-) axis. A 
biaxial anisotropy energy density can be formulated as 
e K = -KM 2 Z + K'Ml where K and K' describe the easy 
and hard axis anisotropy pointing along the z-axis and 
the x-axis, respectively. The case K' — describes a uni- 
axial anisotropy along the wire axis. The above energy 
density is given in units of fioM s 2 , where /zo is the vacuum 
permeability, and M s = \M\ such that the local magne- 
tization M enters here as a unit vector, fh — M/M s . 

If a current is driven through the nanowire, the spatio- 
temporal magnetization dynamics is governed by the so- 
called modified LLG equation with additional AST and 
NAST terms 0, 

dm 



m = — Kim x Ht + am x — (u- V)m+/3m x Uu- V)m], 

dt 

(1) 

where I7I, a, and /3 are the gyromagnetic ratio, the 
Gilbert damping coefficient, and a dimensionless coef- 
ficient describing the NAST strength, respectively. j3 
is usually of the same order as the damping in ferro- 
magnetic metals. The velocity u points along the flow 
direction of the itinerant electrons which is usually the 
wire axis although perpendicular current injection has 
also been proposed [L3|. Thus, V = d/dz in ID, and 
u = g/j,BJP/{2eM s ) where (among standard notation) J 
and P are the density and the spin polarization of the 
current, respectively. 
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The (normalized) effective field h t = H t /M s is given 
by the variational derivative of the total energy density 
(per unit section-area) E = J_ dze(z) with respect to 

magnetization, h t — —SE/Sm(z). The local energy den- 
sity is given byQ 

e(z) = -Km 2 z + K'm 2 x + A(9' 2 + sin 2 9<j)' 2 ) -rh-h, (2) 

where A describes the exchange interaction, and h 
is the normalized external field. Moreover, we have 
adopted the usual spherical coordinates, m(z,i) — 
(sin 9 cos 4>, sin 9 sin (f>, cos 6) where the polar angle 9(z, t) 
and the azimuthal angle 4>{z,t) depend on position and 
time, and the prime denotes spatial differentiation. 

Following Ref. 0, @] , we will focus on DW structures 
fulfilling <p' = (f>" = 0, i.e. all the magnetic moments 
synchronously rotate around the easy axis in space. Then 
the dynamical LLG equations take the form 

T9 = a{2A9" - K sin 29 - K' sin 29 cos 2 (f>) 
+ K' s'm9sm2(j)+ (ah e + h^) - (1 + a(3)u9' , 

T sin 94) = - (2A9" - K sin 29 - K' sin 29 cos 2 </>) 
+ aK' s\i\9si\i2(j)+ (ah^ - h e ) - (a - j3)u9' , (3) 

where we have defined Y = 1 + a 2 and introduced a 
dimensionless time via t n> i|7|jVT s . hi(i = r, 9, </>) are the 
components of the external field in spherical coordinates. 

The linear DW motion under field or ST in a biaxial 
wire has already been discussed in the literature 0, HI- 
Let us first review the case of a flat wire with biaxial 
anisotropy. Following the pioneering work by Schryer and 
Walker [5| we concentrate on solutions fulfilling 4>{z,t) = 
</>o = constant. This assumption implies all the magnetic 
moments move in a plane and is valid at sufficiently low 
fields [5|. Substituting the travelling- wave ansatz tan | = 
exp into Eqs. ([3]), we obtain 

Tv = A(ah - K' sin 20 o ) + (1 + aP)u , 
A(aK' sin 20o + h) - (a - j3)u = , 

where A = y/A/(K + K' cos 2 <p ) is the DW width, and 
the external field with a magnitude h along the z axis 
is assumed (he = —hsm9, = 0). Thus, the constant 
plane angle 0o and the DW velocity satisfy 

aK'shx2<f>o = (a-P)u/A-h, v = Ah/a + /3u/a. (4) 

This solution just recovers the Schryer- Walker result [H 
in the presence of ST and implies /? ^ i.e. NAST is 
indispensable for nonzero DW velocity 0. In particular, 
for | (a — P)u/ A — h\ > aK' the sine in Eq. @ becomes 
larger than unity, and the solution breaks down. There- 
fore, the limit of uniaxial anisotropy K' — > cannot be 
reached within the above travelling- wave ansatz. 

However, the results (j4j, devise the following scheme 
to experimentally determine the spin polarization of 



the current combining measurements of field-driven and 
current-driven DW dynamics: First perform velocity 
measurement using field- and current-driven DWs sep- 
arately and obtain the quantities A min /a = C\ and 
pP/a = C 2 . Here A min = y/A/(K + K') is the min- 
imum DW width and C\ can be obtained by extrap- 
olating the data to h — > 0. Then apply a fixed field 
such that the Walker limit is reached and the DW width 
reaches its maximum, A max = y 1 'A/ (K + K'/2). By in- 
jecting a spin-polarized current and subsequently lower- 
ing the current density and monitoring the decrease of 
DW width one can reach the situation (a — j3)u = Ah 
implying sin 2(po = and again A = A m j„ . (If j3 > a one 
may reverse the injected current direction.) Now using 
(a — f3)P / A m i n = C3 one can infer the current polariza- 
tion P = C1C3 + C 2 and /3 /a = C 2 /(C 1 C 3 + C 2 ). A large 
anisotropy i.e. K' ^ K and the resolution for observ- 
ing the DW width variation are the key points for this 
scheme. 

Now we look at the DW motion in a uniaxial wire 
(K' = 0) which is essentially different from the biax- 
ial case. Using the travelling- wave ansatz tan | = 
exp ( z ~^ t ) where now A = y/A/K, and also for a static 
field h applied along the z axis, one can straightforwardly 
obtain the following expressions for the velocity and the 
precession frequency from Eq.Q, 

Tv = aAh + (1 + a/3)u, T<]> = h-(a- (3)u/A. (5) 

From Eq. |5]), we immediately conclude an important 
result that sustained DW motion is possible in a cylin- 
drical nanowire even in the absence of NAST, a finding 
contradicting leading opinions in the recent literature. 
Moreover, Eq. ([5]) also shows that, at small damping, 
the velocity of a current-driven DW is essentially propor- 
tional to u while in the field-driven case we have v ~ aAh 
and the DW motion is suppressed. Thus, in a uniaxial 
nanowire, a current is more efficient in driving a DW 
than an external field. Furthermore, the precessional fre- 
quency also occurs in the recentl y s tudied DW driven 
electromotive force V em f = ±h<j)/e [14j . which should also 
motivate future experiments. In particular, in a cylindri- 
cal nanowire the current polarization can be easily de- 
termined experimentally by measuring the velocity of a 
current-driven DW. 

Let us now discuss the energy variation in the field 
and/or current driven case in both types of nanowires. In 
general, the change rate of total magnetic energy (from 
anisotropy, exchange interaction and external magnetic 
field) can be expressed as 

/oo 
dzrh ■ [am + u(m x d z rh) + (3ud z rh] 
-00 

= P a + Past + Pn ast ■ (6) 

Here P a = —2aA(4> 2 + v 2 / A 2 ) is the total dissipa- 
tion power due to damping, and Past = —2u(p and 
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-Pnast = 2/3uv I A are the energy pumping rates induced 
by AST and NAST, respectively. Note that the DW 
width A follows different expressions in the biaxial and 
uniaxial case, and v and <p are given by Eqs. (TJ| and ([5]). 
Remarkably, the AST does not contribute to any change 
in energy for a biaxial wire (<f> = 0) below the Walker 
limit. Moreover, it is straightforward to verify that in 
both cases it holds E — —2hv which is just the released 
Zeeman energy per time of a rigid DW traveling with 
velocity v in the external field /i[6(. Thus, the relation 
E = —2hv can also be used as a definition of the DW 
velocity from an energetic point of view. In particular, 
a moving DW driven purely by current (zero magnetic 
field) conserves its energy. 

Let us now consider the change of angular momen- 
tum L = J dzm{z, t) (per section-area and in units of 
— fXoM s /\j\) due to ST. In one hand, the angular momen- 
tum change due to rigid DW motion is L — Tt<j)Ae l f ) + 2vz. 
On the other hand, LLG equation ([1]) gives 

/oo 
dz[h t x rh + am x rh — ud z rh — (3urh x d z rh] 
-oo 

= T pre + T a + Vast + Tnast • (7) 

For our uniaxial wire, T pre — ixhAe^, T a = —ixcme^ + 
2a(j)Az, Tast = 2uz and Tnast = ^(iue^. Hence, the 
AST pumps a longitudinal spin to the wire to push the 
DW propagating while the NAST and external field pro- 
vide the transverse torques forcing DW precession around 
the wire axis. Furthermore, T a due to the DW preces- 
sion (</>) provides an extra longitudinal torque and that 
due to the propagation (v) results in an effective trans- 
verse force. These two damping effects are reminiscent 
of Barnett effect [IB] and Einstein-de Haas effect re- 
spectively. In the present case, both effects originate 
from non-zero damping a ^ 0. Moreover, two of us 
have recently studied DW motion driven by a circularly 
polarized microwave, which also embodies the Barnett 
effect UJ. We also note that Eqs. ([5]) can equivalently 
expressed as 

v = u + a<i)A, <j) = h + (flu - av) / A, (8) 

which shows explicitly that the DW precession con- 
tributes to its velocity as a<j)A, and the DW translation 
contributes to the precession as —av/ A. The above con- 
siderations refer to a uniaxial wire. In the biaxial case, 
however, T pre acquires an additional term from the trans- 
verse anisotropy K' and the DW precession is absent. As 
a result, the effect of AST is canceled and only NAST af- 
fects the DW velocity. 

Before ending the letter, we would like to give a 
practical example of the DW motion in a permal- 
loy uniaxial or biaxial nanowire. We employ the 
standard data M s = 800kA/rrL A = 20nm, a = 
0.02, £ = 0.04, P = 0.4[| fl. For a uniaxial 



wire, w(m/s)» 0.7£f (100Oe)+2.9 J(10 7 A/cm 2 ), </>(GHz) 
1.76iJ(100Oe)+0.003J(10 7 A/cm 2 ); for a flat wire, 
w(m/s)» 1760#(100Oe)+5.8J(10 7 A/cm 2 ) below the 
Walker limit. Thus the current- induced DW velocity 
has the same order in both types of wires. The field- 
driven azimuthal processional frequency in a cylindrical 
wire is in GHz regime while current-driven precessional 
frequency is in MHz regime. 

In summary, we have demonstrated, in the framework 
of the modified Landau-Lifshitz-Gilbert equation, that 
adiabatic spin torque only is enough for current-driven 
DW motion in a cylindrical (uniaxial) nanowire. This 
finding contradicts leading opinions expressed in the re- 
cent literature. Moreover, we have also discussed the mo- 
tion of a rigid DW being subject to (adiabatic or nona- 
diabatic) spin torque in flat or cylindrical wires from an 
energetic and angular momentum point of view. Finally, 
we have also proposed an experimental scheme to mea- 
sure the spin current polarization by combining both field 
and current driven DW motion in a flat (biaxial) wire. 
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